Universal alternating order around impurities in antiferromagnets 
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The study of impurities in antiferromagnets is of considerable interest in condensed matter physics. 
In this paper we address the elementary question of the effect of vacancies on the orientation of 
the surrounding magnetic moments in an antiferromagnet. In the presence of a magnetic field, 
alternating magnetic moments are induced, which can be described by a universal expression that 
is valid in any ordered antiferromagnet and turns out to be independent of temperature over a large 
range. The universality is not destroyed by quantum fluctuation, which is demonstrated by quantum 
Monte Carlo simulations in the two-dimensional Heisenberg antiferromagnet. Physical predictions 
for finite doping are made, which are relevant for experiments probing Knight shifts and the order 
parameter. 

PACS numbers: 75.10.Jm, 74.25. Nf, 75.20.Hr, 75.40.Mg 
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The intentional doping of antiferromagnetic materials ^ ^ .r-^.^^.^..-- 
has become a useful tool in order to study the compli- — 
cated physics in the context of high temperature super- ^ % ^ 

conductivity and quantum magnetism 0, S, 0j S Si- ^ ^ 

Large alternating magnetic moments around static non- ^ ^ ^ -..^ ^ 

magnetic impurities are observed in Kiiight shift experi- 

ments when a uniform field is appliediiii. Theoret- ^, (Color online) Schematic illustration of the effect of a 

ical studies have shown that vacancies m low-dimensional ^^^^^^^ ^^^^^^^ antiferromagnet. The spins "cant" with 

an angle & towards the field corresponding to a small uniform 
magnetization. Due to the broken sub-lattice symmetry the 
order may be "tilted" by an angle a relative to the plane 
normal to the field corresponding to an induced alternating 
magnetization around the impurity. 



antiferromagnetic backgrounds give rise to locally en- 



rise to locally en- 
hanced antiferromagnetic correlations 0, [E B El, 
13l . Il4j , which strongly depend on the microscopic model 
and temperature in the low dimensional models. 

In this work, we show that in generic ordered antifer- 
romagnets the alternating local moments in the vicinity 
of vacancies can be quantitatively described by a univer- 
sal expression which only depends on the field B, but is 
surprisingly independent of temperature, quantum fluc- 
tuations, and microscopic details. The mechanism which 
gives rise to the alternating moments is a local tilting of 
the order parameter due to the broken sub-lattice sym- 
metry by impurities. In contrast to the pure sample, 
where the order parameter is always confined in the plane 
normal to the field, a large alternating order parallel to 
the field is induced as schematically depicted in Fig. [T] 
The calculations agree remarkably well with quantum 
Monte Carlo (QMC) simulations without any adjustable 
parameters even in two dimensions D = 2, where quan- 
tum fluctuations are strongest. 

The typical antiferromagnetic Hamiltonian 



(1) 



describes the magnetic behavior realistically even for 
rather complex materials despite its simplicity. We con- 
sider systems with bipartite lattices of dimension D >2, 
where the sum in Eq. ([1]) runs over nearest neighbor sites. 
Generically, the dominant interaction J > comes from 
the Coulomb forces via the exchange mechanism and is 



therefore isotropic. The rotational symmetry is broken 
by applied and crystal fields B in units of g/is, which 
are typically small compared to the interaction B < J. 
The direction of the field defines the z-axis of our coordi- 
nate system, which does not need to coincide with any of 
the lattice directions. For bipartite lattices of dimension 
D > 2 the model system ([1]) is known to have finite Neel 
order at sufficiently low temperatures for both quantum 
and classical spins Si of any size s[3l. The ordered state 
remains stable over a large range of perturbations by im- 
purities and frustrating interactions. 

In order to obtain an intuitive picture of the physical 
behavior, let us first consider a highly simplifled model of 
the Hamiltonian ([T]). The long-range order spreads over 
the entire sample, so it might be justified to describe all 
ordered spins on one sub-lattice A by a common direc- 
tion fiA — {sill 9 A sin 0^, sin 0^ cos0yi, cos6'^) and anal- 
ogously for sub-lattice B. In this case, the interaction 
energy is always minimized by a relative azimuthal angle 
(t^A — 4>B = TT, SO that the effective energy is given just in 
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terms of the polar angles 

Ecs = JzNs'^fiA ■ fiB - sNB{n''j^ + n%) 

= JzNs'^ cos{9a + 9b)- sNB{cos 9a+cos9b) 
= JzNs^{2sm^ 5 - I) - 2sNBcosasuv6 (2) 

where z is the number of nearest neighbors and N is 
the total number of sub-lattice sites in the sample. The 
angle 5 = {tt — 9a — 9b) /2 corresponds to a uniform 
"canting" of all spins on both sub-lattices towards an 
applied magnetic field as shown in Fig.[TJ The angle a — 
{9b — 9a) measures the "tilt" of the antiferromagnetic 
order relative to the plane that is normal to the field. 

Below saturation B < 2szJ the energy is minimized 
by setting sin 6 — i^pf cos a which gives an effective low 
energy description for a 
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NB^ 
2zJ 



■ sm a 



En 



(3) 



where Eq = - JNs^z - NB'^/2zJ. The physical inter- 
pretation of this simple model is textbook knowledge [l6|: 
All spins align slightly towards the magnetic field m = 
s(sin5) = Bxa. cos a with a susceptibility xj. ~ l/2zJ 
that is largest when the magnetic field is perpendicular 
to the order and therefore there is a small energy gain for 
the Neel order to be in the plane normal to B (i.e. a = 0). 
Since the energy gain is small, the order may point in an- 
other direction in realistic materials where the sub-lattice 
symmetry is broken. A common source of sub-lattice 
symmetry breaking is disorder and impurities which is 
the topic of this paper. 

Let us first consider a single vacancy in the frame- 
work of the simple model above by reducing the size of 
the corresponding sub-lattice vector by one spin = 
s{N — \)hA- Starting from Eq. ^ the susceptibility for 
5 remains the same for large N . However, the size of the 
two sub-lattice spins is not equally large and therefore a 
net coupling to the field remains in the effective energy 
as a function of the z-component of the alternating order 



sma 



E^s{n,) - NB^xi-nl ~ sBn, + Eq 



(4) 



where we have also used that the dependence on S is small 
and irrelevant in the direct coupling term. Even though 
N is large, the second term will ensure that the expec- 
tation value of the impurity induced alternating order 
along the field is always non-zero 



BNx± \2 j\e—^B^«x±Pdx^ 



(5) 



where we have assumed the thermodynamic limit N ^ 
fiJ {p = l/fc^T). In the limit of large and small fields. 



respectively, we find 
sB/3T 
S/2NBXJ 



for Nx±B^ < T 
for Nx±B^ > T 



(6) 



In the first case of very small fields, the alternating 
response to a uniform field is described by a classical 
susceptibility, which also directly follows from Eq. (|4]) if 
only the second term is kept (i.e. x±B 0). Therefore 
a tilting of the order parameter of order a ~ B/T is ex- 
pected which is larger than the uniform canting 5 ^ B / J 
in the ordered phase T <^ J. By QMC simulations it 
was shown that a corresponding alternating order is in- 
duced throughout the lattice by a single vacancy in the 
limit of linear response llj , which is consistent with the 
assumption that a describes the tilting of all spins. The 
corresponding impurity susceptibility is given by a clas- 
sical Curie behavior s^/ST, as first predicted in Rcf. ^ 
and confirmed by QMC simulations in Ref. JJ} in the 
limit of linear response Nx±B^ <C T. This limit is only 
relevant in the case where the domain size N is restricted 
by disorder or boundaries. 

However, if N is macroscopic, the limit Nx±B^ ^ T is 
already reached for any naturally occurring background 
field, so that the second case in Eq. ^ is the more in- 
teresting limit for the description of realistic impurity 
effects. The induced alternating magnetization decreases 
with increasing field and the behavior is independent of 
temperature since corrections from the second term in 
Eq. ([5]) are exponentially small in the macroscopically 
large scaling variable Nx±B^f3. This remarkable behav- 
ior will survive even in more refined models and give rise 
to a universal temperature independent description as we 
will see. 

In order to make the model more realistic, the angle 
a can be interpreted as a local tilting that is dependent 
on position in order to reflect the fact that the first term 
in Eq. ^ is an effective potential that acts on all spins 
in the lattice, while the second term arises from the va- 
cancy locally at the origin. There is an energy cost to 
change the direction of the order parameter from one 
lattice site to a neighboring lattice site corresponding to 
the so-called spin-stiffness ps, so that Eq. ^ has to be 
generalized to an energy functional for 



i?K(r)] 



d^r (|(Vn4r))V^S^n^(r)) 
sBn,{0), (7) 



where we have replaced the sum over both sub-lattices 
by an integral for convenience. The energy density in 
the first term is reminiscent of the non-linear sigma 
model 18, U^, but only for one component and with- 
out the imaginary time direction describing the quantum 
fluctuations. 
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In order to calculate the expectation value of {nz{ro)) 
at any position ro , it is useful to define a generating par- 
tition function 



= / PK(r)] exp i-PE [n,{r)] + 7n.(ro)) • (8) 



The expectation value is then given by the logarithmic 
derivative 

(f^^(ro)) = dylnZ^\^=o. (9) 
In momentum space the partition function becomes 
Z.,^ /■p[„^(q)]e[/'^"q("/3£,|«4q)l^+/,(7)«4q))] (10) 



where 



(11) 



/q(7) = (/3Bs + 7COsq-ro)/(27r)^/2. (12) 
The expectation value is therefore 
(nziro)) = d^lnZ^l^^Q 



rf^q i3scos(q-ro) 
(27r)^ ipsq^+X±B^) 

s E 
2-iTp 



7=0 



s B ^-Br/c 



D^2 

D = 3 



(13) 
(14) 



where c = \fpslxL is known as the spin-wave veloc- 
ity and Kq is the modified Bessel function. This result 
is remarkable in two ways: first of all it turns out to 
be completely independent of temperature and secondly 
it is independent of the underlying detailed microscopic 
model. Therefore, the formula in Eq. can be taken 
as a universal description for all antiferromagnets in the 
ordered phase. Variations in the lattice structure, frus- 
tration, quantum effects, and the detailed microscopic 
parameters only renormalize the spin stiffness and the 
uniform susceptibility X-Lj but not the functional behav- 
ior in Eq. (HH). For spins close to the vacancy r > 1 the 
tilting Uz ^ sB /A'Kps remains typically less than satura- 
tion, but larger than the uniform canting a > 5, so that 
spins on the same sublattice as the vacancy tend to align 
against the field. 

It can be checked that the functions in Eq. are 
solutions of the diffusion equation 



(15) 



that also follows from minimizing the energy functional 
Eq. ([7]). In lattices where the spin-stiffness is not 



isotropic, the result can be generalized by taking ps as 
an anisotropic diffusion coefficient. 

As a concrete example, we will now consider the spin- 
1/2 Heisenberg model on a 2D square lattice, which is 
possibly the most studied antiferromagnetic model, since 
it has received much attention in connection with high 
temperature superconductivity, but also because it is 
an interesting case where quantum fluctuations strongly 
compete with Neel order. 

In Monte Carlo simulations we have used the stochastic 
series expansion with directed loop updates [lO] in order 
to extract the magnetic moments {iriz) around a single 
vacancy in small magnetic fields in the ordered phase 
^(T) ^ L (here L = 128) as shown in the inner inset of 
Fig. [21 In the plane perpendicular to the field the order is 
fluctuating, so that {rrix) — {my) = 0. For the moments 
parallel to the field, we expect to find a large staggered 
magnetization according to Eq. pi)) 



(m,(r)) = {-lym 



max 



"27rps 



B 



(16) 



in addition to the less interesting small uniform canting 
5. Here Wmax ~ 0.308 is the maximum order in the 2D 
Heisenberg model which is reduced from s = 1/2 due to 
quantum fluctuations. In fact, there are no adjustable 
parameters in Eq. p6p since all relevant parameters have 
long been established to high precision by independent 
methods [2l|: 



XA 



« 0.308, Ps ~ 0.18J, 

0.065/ J, c = Vp./Xi ~ 1-67J. (17) 



By extrapolating the numerical data for ruz (r) on the 
even and the odd sub-lattice separately and taking half 
the difference, we extracted the staggered magnetization 
rriait (r) around a vacancy. The resulting alternating am- 
plitude TOait is completely isotropic and can be plotted 
as a function of the geometrical distance r — |r| only 
as shown in Fig. [5] for different fields and temperatures. 
While the size of mait is proportional to the field, the 
drop-off is shortened for higher fields so that the inte- 
grated amplitude decreases with increasing field as also 
reflected in the simple model of Eq. ^ . The agreement 
with Eq. (|16p is remarkably good even on a logarithmic 
scale and for widely different fields and temperatures, 
which we take as confirmation for the general validity of 
the result in Eq. (fT3|l . Since there were no adjustable 
parameters, we conclude that the quantitative predictive 
power for static expectation values of the hydrodynamic 
model in Eq. ([7]) is not changed by quantum fluctua- 
tions. From a theoretical point of view this means that 
the renormalized classical model[3, [3| can be taken for 
quantitative calculations anywhere in the antiferromag- 
netic phase, while microscopic details only affect the val- 
ues of the constants in Eq. p7p . In particular, close to 
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20 25 
geometrical distance i- 

FIG. 2: (Color online) The alternating response mait as a 
function of geometrical distance r from a vacancy in a 2DHAF 
at different fields and T — 0.025J from QMC simulations 
compared to the universal theoretical prediction mmax^z in 
Eq. ((16} without any adjustable parameters (dashed black 
lines). Inset: mait at _B = 0.2J. Even at higher T — O.IJ 
no deviations from Eq. (|16p can be seen on a logarithmic 
scale. At still higher T the induced order first increases for 
0-3J > T > 0.2,1 and then decreases again for T > 0.3J. 
Inner inset: Magnetic moments at B = O.IJ and T — O.IJ 
alternating between the extrapolated amplitudes on even and 
odd sublattices. 



a critical point ps and TOmax become vanishingly small, 
but the model remains valid. 

A breakdown of the universal formula in Eq. (|13p must 
occur at the transition temperature to the disordered 
phase. In the 2D simulations we find indeed that any 
temperature dependence is exponentially supressed un- 
til the Kosterlitz-Thouless temperature is approached 22 
Tkt ^ 0.2J. However, the induced alternating magneti- 
zation is surprisingly enhanced by increasing temperature 
near Tkt as shown in the inset of Fig. [21 Only at still 
higher temperatures T > 0.3 J the induced order is finally 
reduced as expected, leading to a non-monotonic behav- 
ior with temperature. While we have no explanation of 
this exotic effect in terms of our model, we hope that 
future works on this topic may uncover this mystery. 

Finally, we would like to generalize our results to finite 
impurity concentrations p. For higher fields/small con- 
centrations p < {B/c)^ the impurities are sufficiently far 
apart to be treated independently (dilute limit). In this 
case, the above conclusions are unchanged and the mag- 
netic order is tilted locally in the vicinity of each vacancy. 
At smaller fields/larger concentrations p > [B/c)^ the 
impurities become correlated and enhance/annihilate the 
tilting effect depending if they are on the same/opposite 
sub-lattices ll| . In this disorder limit all impurities be- 
come correlated and the tilting is again nearly uniform 
throughout the lattice, with an effective total impurity 
strength that is given by the difference of the vacant 



sites on each sub-lattice \Na — Nb\ ~ ^/ pN in a do- 
main of size N . In this case the simple model in Eq. (j4|) 
remains valid with the effective size of the spin s in 
the last term replaced by sy/pN. The average univer- 
sal tilt in Eq. ^ throughout the domain is then given 
by m 

max'-' 

In conclusion, we have analyzed the induced alternat- 
ing magnetization around vacancies in ordered antifer- 
romagnets in quantitative detail. Large alternating mo- 
ments are induced parallel to the field, which corresponds 
to a tilting of the order parameter. The induced order 
decays with distance at a rate that is independent of 
temperature and inversely proportional to the field c/ B. 
From a theoretical point of view we have demonstrated 
that the renormalized classical description gives an in- 
tuitive insight into the mechanism on how the alternat- 
ing magnetization arises. At the same time, the the- 
ory gives good quantitative agreement with QMC simu- 
lations even in 2D where quantum fluctuations are large. 
At large impurity densities p > {B/c)^ the impurities 
become correlated and give rise to a tilt of the order 
throughout the sample towards the field direction. The 
predicted effects can be observed by Knight shift exper- 
iments like NMR and /iSR, or by investigating the or- 
der directly via magnetic neutron scattering. Numeri- 
cally we find an enhancement of the induced order near 
the Kosterlitz-Thouless transition Tkt ^ 0.2J which is 
counter-intuitive and calls for further investigation. 
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ration was supported by the Nordforsk Network on Low 
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